Thermodynamics of strongly frustrated magnet in a field: Ising antiferromagnet on 

triangular Husimi lattice 
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Some strongly frustrated magnets such as the "spin-ice" compounds fail to produce any magnetic 
order at finite temperatures even in the presence of magnetic field. Still they have very unusual low-T 
thermodynamic properties related to the field-induced ground state transitions. Here we show that 
general qualitative picture of such peculiar thermodynamics can be obtained in the antiferromagnetic 
Ising model on the triangular Husimi lattice. The analytical results for this model show magnetic 
plateaus, entropy spikes, crossing points and peculiarities in magnetic susceptibility and specific heat 
behavior reflecting the existence of ground state transitions. These signatures of strong frustration 
may help in search of new frustrated magnets and in the interpretation of experimental data. 

PACS numbers: 75.10.-b, 75.40.Cx 
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The magnetic ions in some crystal lattices do not have 
ordered magnetic phases at finite temperatures in spite 
of magnetic interactions between them. The origin of 
this phenomenon hes in the specific topological proper- 
ties of such lattices which forbid the existence of unique 
ground state configuration of local spins. Instead the 
spin systems on these "geometrically frustrated" lattices 
have numerous ground states the number of which grows 
exponentially with system size [l|-|3|. Some symmetry 
breaking perturbations (magnetic field, anisotropic stress 
etc.) may lift this degeneracy and give rise to some exotic 
ordered or partially ordered phases [4J, |5| . 

Yet in some strongly frustrated magnets such per- 
turbations fail to produce any magnetic order at finite 
T. The notorious example is the "spin-ice" compounds 
with pyrochlore lattices consisting of the corner-shearing 
tetrahedra [6|. Such magnetic systems may only have 
zero-temperature transitions under magnetic field vari- 
ations between some degenerate ground states. Some 
signatures of these transitions can be found in the low- 
temperature thermodynamics. Thus magnetic plateaus, 
entropy spikes and crossing points shows up in the model 
of spin-ice compounds ^ . Yet the theoretical description 
of such effects on realistic 3d lattices needs considerable 
efforts even in cases of special field orientations where the 
model can be mapped to the frustrated Ising antiferro- 
magnet [7|. 

So it may be useful to consider simplified models to 
gain the general insight into the origin and overall pic- 
ture of the peculiar low-temperature thermodynamics of 
strongly frustrated magnets. Here we show that the 
model of frustrated Ising antiferromagnet on triangular 
Husimi lattice can qualitatively reproduce all peculiari- 
ties of spin-ice model in [111] field [7| which are related 
to the field-induced ground state phase transitions. This 
is quite simple model allowing for the analytical descrip- 
tion of all thermodynamic parameters thus giving the 
consistent picture of the frustration-induced anomalies. 



So these results may help in the interpretation of exper- 
imental data and in search for new strongly frustrated 
magnets. 

Fragment of infinite triangle Husimi lattice composed 
of Husimi trees is shown in Fig. [T] We consider Ising anti- 
ferromagnet on this lattice. Its partition function can be 
obtained through consideration of the recursion relations 
for partial partition function of Husimi tree summed over 
all spins except the root one 
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FIG. 1: Fragments of Husimi lattice (a) and Husimi tree(b). 
Infinite species are obtaned by sequential addition of triangles 
to outer sites. 
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Here Tr denotes the summation over spins si ands2, K 
J IT, h = H/T. Defining the effective field 



2 Z„(-l) 



we get for it from Eq. [T] 



^1 cosh2{fn + h-K) + e^^ 
"^"+' ~ 2 "" cosh 2 (/„ + h + K)-^e^^ 



(2) 



The stationary point of Eq. [2] foo — f {K, h) defines 
all thermodynamics of the model. Indeed, with foo = 
f [K, h) we can obtain the free energy per spin F (T, H) 



as follows. Let us consider 12 Husimi trees and divide 
them in 4 triples. Adding to each triple 3 bonds we can 
get 4 Husimi lattices. Otherwise we can divide 12 trees 
in 3 quadruples and then adding to each quadruple 1 site 
and 6 bonds we get 3 Husimi lattices , cf. Fig. [5] In 
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In the low temperature region Eqs. (6), (7) can be fur- 
ther simplified. First, when field H > is also small, 
max {H, T) < J, (y > z > 1) we have 

zD^ + 2D^ - 2zD -1 = (8) 

3F/T = -2K + 3 ln(;z) + In (l + D^) - 

2 In (z^D^ + 2zi:i) . (9) 

Thus D and F/T in this region depends only on z and 
this results in very unusual behavior of thermodynamic 
variables when H,T -^ 0. Thus magnetization ([5)), en- 
tropy 



FIG. 2: Construction of Husimi lattice from 4 Husimi trees. 

the infinite lattice limit the free energies of all Husimi 
lattices in both constructions will tend to the similar ex- 
tensive values with the only difference that in the latter 
construction we have 3 more sites than in the former one. 
Thus for the free energy per spin F (T, H) we get 

3i^ = -3T lim {ln[/„}+4T lim {InK} 



Un^Tr Y[ Zn {si)e 

3 

Vn = Tr n Zn {Si)e 

1=1 



{h — KsQ)(si+S2+S3 + S4) + hso-K{siS2+S3S4,) 
h{si+S2+S3)-K{siS2+S3Si+S3S2) 



Performing here the summation over spin values and 
using the stationary point equation for foo = f (K, h) 

1 cosh2[f + h-K) + e'^ 
■' 2 "^ cosh2 {f + h + K) + e^^ ^ ' 

and definition of / we finally get 

3F/T = 2K~\n2-2f + \n cosh {h + 2/) - 

21n[cosh2(/i + / + X) + e2^] (4) 

Differentiating F with respect to H and T we can get 
magnetization, magnetic susceptibility, entropy and spe- 
cific heat of the model. In particular, we get the simple 
expression for magnetization 

m = -—= tanh (2/ + h) = ;^^^, D = exp (2/ + h) . 

(5) 
Introduction of new variable D instead of / is useful be- 
cause Eq. 3 becomes with it just simple cubic equation 



yzD'+{2y-y-'z^)D^ + {y-'z-' 



2yz)D-y = (6) 



Here y = exp 2K, z = exp h. Note that as m > 
at _ff > physically relevant solution for D should be 
greater than 1 for positive fields which we consider fur- 
ther. Expressed in these new variables free energy (|3]) 
has the form 

3F/T ^3\n{zy) + \n (l + D^) -2\n{z^y^D^ + 2zy^D + l) 

(7) 



S = -—^-ln (1 + 2^^)' 
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reduced magnetic susceptibility 
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(10) 
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and specific heat 



C ^T 



dS_ 
df 



(Inz)^x' 



(12) 



are the functions of the ratio H/T only. So they are 
constant along the lines H/T = const, when iJ, T — ^ 0. 

The physical solution {D > 1) to Eq. [5] (which is also 
the only stable one with x > 0) h^s the form 



D{z) 



3z 



Re 



--z^-8 + i-zR(z) 
2 2 ^ ^ 



R{z) = V3(32z4 + 61z2 + 32) 

Here the power -1 means the main branch of the power 
function. This D(z) grows monotonously from D{1) = 1 
to D (oo) = v^- So at if = we have 

TO = 0, x = 4/5r, 5*= iln^ «0.501, C = 

while at T = 

m = l/3, x = 0, 5 = iln^«0.096, C = 0. 

In general these thermodynamic variables has different 
limiting values at 77, T — >■ depending on the way to 
this point in the _ff , T plane. The general picture of their 
behavior near i/ = T = is shown in Fig. [3] 

Similar scaling features can be found in strong fields 
H > 2J > T (z > y > 1). Here we have 



D = -{x+ \Jx'^ + \ 



X = zy — exp (h — AK) 



3F/T= -2K + lnx + ln{l + D^) -AlnD 
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FIG. 3: (color online) Field and temperature dependences of 
thermodynamic parameters at H,T <^ J. 



FIG. 4: (color online) Field and temperature dependences of 
thermodynamic parameters at T <C J, H ~ 4J. 



So D and F/T in this region also depends only on one 
variable x (or {H -- AJ) /T). Accordingly, this scaling 
holds for the thermodynamic quantities 
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C= (Inx)^x' 



Their behavior near iJ == 4 J is shown in Fig. SI At T = 
we get 



1 
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i7<4J, m=^, S=^\tl^, x' = C' = 0, 

ff = 4J, m = 2 5 = 1 In M w 0.388, x' = #, C = 

' 5 ' 3 5 ' -^ 75 ' 

H>U, TO = 1, 5 = x' = C* = 

Note that the results for T = 0, if < 4J coincide with 
those for T = 0, if < J. So at < if < 4 J we have a 
plateau in the field dependence of magnetization [m — \) 
and other thermodynamic parameters and at if > 4J 
there is another plateau. They are shown in Fig. [5] for 
TO and 5. 

Thus at T = we have two field-induced first-order 
transitions at if = and at if = 4 J. The nature of these 
transitions is quite apparent. In zero field we have highly 
degenerate ground states having two parallel spins and 
the anti-parallel one in each triangle. Their number can 
be computed using Pauling- Anderson- type estimates [8j . 
First we consider the spins as belonging to A^/3 indepen- 
dent triangles each having 6 ground state configurations. 
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FIG. 5: (color online) Field dependences of magnetization 
(solid lines) and entropy (dotted lines) at T = and T = O.IJ. 



Each site in such configurations takes the values -1-1 and 
-1 with probability 2- Hence the probability p that the 
bond triangles connecting the independent ones have also 
the lowest energy configurations is p = '^ (5) • Thus 

the number of ground states is F = (6p) — (|) 
giving the quoted above result for the zero-field entropy 



S = N~^ InF. Thus Pauling- Anderson entropy estimate 
[8| neglecting the correlations between spin ordering in 
the next nearest triangles is exact for zero-field ground 
states on the Husimi lattice. 

The arbitrary small field lifts partially this degeneracy. 
Here we can also apply the Pauling- Anderson approach. 
Now only three triangle configurations (permutations of 
(1, 1,-1)) give the lowest energy. In such configurations 
each spin is 1 with probability 2/3 and -1 with proba- 
bility 1/3. So the probability p that the bond triangles 
connecting the independent ones also have ground state 
configurations is now p = 3i(|) =| and the number 

of ground states is F = (3p) —(f) i^^ accordance 
with above S ioi < H < 4 J. So here we have another 
phase with lower S and m ~ ^. 

One more phase (ferromagnetic) appears at iJ > 4J 
with 5 = and m — 1 while at iJ = 4J special 
ground states exist. At this point the energy of ferromag- 
netic configuration of two triangles with one common site 
(dashed lines in Fig. [2]) does not actually depend on the 
direction of the central spin. The magnetization of these 
states TO = I shows that they have the fraction 1/5 of 
spins pointed opposite to the field direction. So we may 
conclude that here the ground state configurations can be 
obtained by the divisions of lattice into couples of trian- 
gles with common site having -1 spin and other four sites 

with +1 spins. Then at 7? = 4 J F = e^^ = {f)^^^ can 
be identified with the number of divisions of Husimi lat- 
tice into a couples of connected triangles. Yet here we can 
not prove this correspondence. We may only observe that 
the huge degeneracy of ground states at ff = 4 J give rise 
to the entropy spike which transforms into sharp maxi- 
mum at finite T, cf. FigjSl quite similar to that found in 
the model for spin-ice pyrochlores [7|. 

In general, the low-T features of thermodynamic pa- 
rameters in a field stem from the existence of T = phase 
transitions. Thus all magnetization curves for different 
low-T values will have to, = | at iJ = 4J so this is the 
crossing point for these curves. This crossing effect has 
been found in the spin- ice model [Tj. But the present 
model shows that it should also exist for all dimension- 
less parameters S, C and x'^ at all small T ^ J they 
acquire their ground state values at H — 4J. Note also 
that the divergences in temperature dependencies of sus- 
ceptibility X ~ T~^ appear only at the field transitions' 



points H — and H = AJ, otherwise x — >■ at T — >• 
coming over a broad maximum. Meanwhile the magnetic 
specific heat C is zero at H = and H = AJ which give 
rise to its double-peaked field dependencies around these 
points (cf. Figs. [31 S]). Along with the entropy spikes 
and magnetization plateaus this may serve as indication 
of the ground state phase transitions in which ground 
state degeneracy is lifted by an external field. 

To conclude we may state that present model of Ising 
antiferromagnet on Husimi lattice can give qualitatively 
adequate picture of the specific low-temperature thermo- 
dynamics which may exist in real 3d frustrated magnets 
when their spin system allows for the Ising-type descrip- 
tion. Apparantly the model with effective dimension 
d = oo may give only approximate information on the 
structure of the ground state configurations in real crys- 
tals and may grossly overestimate their degeneracy as 
compared with the lattices in 3d Eucledean space. This 
fault can be partially remedied by considering the gen- 
eralizations of Husimi lattice where the basic building 
blocks are the fragments of 3d lattices instead of trian- 
gles. Yet failing to reproduce quantitatively the values 
of thermodynamic variables in real lattices it may cor- 
rectly describe the thermodynamic anomalies near field- 
induced phase transitions because of their first-order na- 
ture allowing for the mean-field description. Also the 
consideration of anisotropic Heisenberg model on such 
lattice may provide useful qualitative information on the 
properties of real geometrically frustrated magnets. 
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